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Abstract 

In  pursuit  of  a  more  efficient  and  effective  fuel-air  mixing  for  a  SCRAMjet  combustor, 
this  study  investigated  relative  near  field  effects  of  molecular  weight  on  mixing  and 
penetration  of  different  gaseous  injection  into  a  supersonic  (M=2.9)  cross  flow.  Hehum  and 
argon  gas  were  chosen  as  injectants  because  of  the  large  differences  in  molecular  weights. 
Also,  miytng  enhancement  was  observed  by  injecting  the  traverse  gas  jet  parallel  to  the 
compression  face  of  a  ramp.  Color  schheren  photography  was  used  to  identify  the  shock 
stmctures  and  interactions  in  the  flow  field.  Measurements  of  mean  flow  properties  were 
used  to  estabhsh  the  jet  plume  size,  penetration,  and  concentration  and  to  quantify  the  total 
pressure  loss.  Results  indicate  greater  mixing  and  plume  expansion  can  be  achieved  with 
hehum  compared  to  argon. 
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EXPERIMENTAL  INVESTIGATION  OF  THE  INFLUENCE  OF  MOLECULAR 
WEIGHT  ON  MIXING  AND  PENETRATION  IN  SUPERSONIC  DISSIMILAR 
GASEOUS  INJECTION  INTO  A  SUPERSONIC  CROSS-FLOW 

1.  Introduction 

1.1  Motivation 

The  early  decades  of  aviation  were  characteri2ed  by  the  pursuit  of  speed  and  altitude.  A 
major  push  for  higher  speed  was  initiated  by  the  X-series  of  aircraft  that  explored  the 
supersonic  flight  regime  in  the  1950s  and  1960s. ^  Since  that  time,  air-breathing  propulsion 
systems  have  pushed  material  and  chemical  properties  to  the  limits.  Because  of  these  limits, 
progress  in  the  field  has  been  evolutionary  rather  than  revolutionary  since  the  early  1960s.^ 
Material  sciences  and  computational  resources  have  allowed  systems  to  become  more 
efficient,  but  the  great  strides  in  the  "higher  and  faster"  progression  seen  in  the  late  1 950s 
have  stopped.  The  SR-71,  a  late  1950s  design,  stiU  holds  the  speed  records  for  air-breathing 
manned  aircraft  at  Mach  3+.  The  unmanned  D-21  was  ramjet  powered  and  achieved  Mach 
5+.  The  experimental  X-15,  while  achieving  nearly  Mach  7,  was  rocket  powered.^  The 
benefits  of  a  hypersonic  air-breathing  aircraft  are  obvious.  By  utilizing  the  oxygen  in  the 
atmosphere,  the  aircraft  only  needs  to  carry  the  fuel,  therefore  increasing  the  payload  or 
range.  The  speed  of  such  aircraft  is  also  desired  by  planners  for  fast  worldwide  coverage  and 
as  a  platform  for  orbit  insertion.  These  benefits  have  been  shaped  in  recent  years  into 
requirements  for  the  National  Aerospace  Plane  and  currently  the  "Military  Spaceplane.""* 

For  speeds  in  excess  of  Mach  4,  a  transition  from  turbofan  or  turbojet  to  the  ramjet 
mode  of  propulsion  is  necessary.  The  characteristics  of  the  conventional  subsonic 
combusting  ramjet  are  well  understood.  However,  the  performance  of  the  subsonic 
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combusting  ramjet  deteriorates  at  hypersonic  speeds  due  to  real  gas  effects  and  internal 
loses,  and  the  structural  integrity  of  the  engine  becomes  very  difficult  or  impossible  to 
maintain  due  to  the  high  internal  temperatures  and  pressures.  Thus,  it  is  necessary  to  utilize 
the  supersonic  combustion  engine  (SCRAMjet)  for  such  high  speeds. 

The  SCRAMjet  engine  potentially  offers  outstanding  specific  impulse  performance  at 
hypersonic  Mach  numbers.  The  loss  of  energy  in  the  flow  through  the  engine  is  less  than  in 
a  ramjet  engine,  but  supersonic  combustion  represents  a  major  technological  challenge  in 
engine  development.  The  existence  of  supersonic  throughflow  presents  the  problem  of  how 
to  manage  efficient  and  effective  mixing  and  burniag  of  fuel  for  the  extremely  short 
residence  time  the  working  fluid  is  in  the  combustion  chamber. 

1.2  Background 

The  short  residence  time  of  flow  inside  the  combustor  engine  presents  a  problem  to 
designers  which  must  be  overcome  by  enhancing  the  mixing  and  penetration  of  the  fuel  in 
the  oxidizer.  To  aide  in  the  mixing  and  penetration  of  an  injected  plume  into  a  supersonic 
freestream,  several  schemes  have  been  studied.  Early  studies  focused  on  the  effects  of 
underexpanded  jets  injected  normally  into  the  flow.^>^  While  laying  the  foundation  of  today’s 
computational  approaches,  these  initial  efforts  provided  fundamental  analyses  of  the 
governing  equations.  The  pressure  differential  across  the  front  to  rear  face  of  the  emerging 
jet  was  found  to  enhance  mixing.^’®  In  addition  to  transverse  injection, more  recent 
studies  have  focused  on  retaining  the  jet  momentum  in  the  direction  of  thrust  through  low 
angle  injection  over  or  behind  ramps.^’^^’^^’^^’^^  These  configurations  were  designed  to 
generate  vorticity  from  the  spillage  of  high  pressure  air  off  and  around  the  leading 
compression  surface.^^ 
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Along  with  the  study  of  mixing  enhancements,  a  parallel  effort  has  been  made  to 
quantify  the  concentrations  of  binary  gas  mixtures.  Many  of  these  studies  have  focused  on 
using  a  hot-wire  probe  to  make  calibration  curves  for  known  concentrations  and 
pressures.^^’^^>^^ 

L3  Scope  of  Present  Study 

The  primary  objective  of  this  study  is  to  compare  the  differences  molecular  weight  has 
on  the  mixing  and  penetration  of  an  injectant  gas  into  a  supersonic  crossflow.  This  is 
accomplished  by  examination  of  the  near  field  mixing  of  Mach  1.87  underexpanded 
supersonic  helium  and  argon  jets  transversely  injected  25®  into  a  nominally  Mach  2.9  air 
freestream.  Also,  a  comparison  is  made  between  a  simple  transverse  injection  and  injection 
over  a  penetration  and  mixing  enhancement  ramp. 

Helium  and  argon  were  chosen  for  the  large  differences  in  molecular  weights  between 
the  two.  Helium  has  a  molecular  weight  of  approximately  4  while  argon  has  a  molecular 
weight  of  approximately  40.  Also,  these  weights  approximate  the  trends  found  in  the  two 
most  likely  fuels,  hydrogen  and  hydrocarbons,  for  SCRAMjet  applications.  An  extremely 
light  fuel  and  a  fuel  with  a  weight  an  order  of  magnitude  greater. 

1.4  Limitations  of  Present  Study 

Limitations  of  this  study  include  availability  of  facilities  and  equipment,  number  of 
downstream  locations  sampled,  and  available  data.  A  time  constraint  imposed  on  the  use  of 
shared  resources  limited  the  mean  flow  and  concentration  measurements  to  one  axial  station 
downstream  of  the  injection  site.  The  limitation  in  downstream  test  location  is  partially 
overcome  by  the  fact  the  fuel  plume  under  proposed  flight  conditions  of  a  hypersonic 
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vehicle  will  most  likely  combust  m  the  near-field  due  to  the  extreme  high  temperatures  to  be 
expected  in  the  flowd^  These  limitations  on  the  data  collected  affect  the  understanding  of 
the  detail  mechanisms  acting  in  the  flow  and  more  analyses  of  the  near-field  is  needed  before 
a  complete  understanding  is  attained. 
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2.  Background 

There  have  been  many  experimental  and  computational  studies  in  the  field  of  transverse 
supersonic  injection  and  gas  sampling.  This  chapter  reviews  some  of  this  research  from  past 
investigations.  The  topics  covered  are  flow  and  injection  fundamentals,  transverse  injection, 
vorticity  generation,  and  gas  sampling. 

2.1  Flow  Fundamentals 

Several  flow  variables  and  injection  parameters  have  been  estabUshed  to  aide  in  the 
comparison  of  flows  at  different  conditions.  The  injection  parameters  include  injector-to- 
freestream  ratios  such  as  the  velocity  ratio,  mass  flux  ratio,  and  expansion  ratio.^  The 
dynamic  pressxxre  ratio,  q  ,  is  another  injector-to-freestteam  ratio  that  is  also  known  as  the 
momentum  flux  ratio.  The  primary  variables  of  concern  for  this  study  are  the  injector 
diameter  based  Reynolds  number,  d,  the  boundary  layer  depth  to  injector  diameter  ratio,  7 , 
and  the  flow  total  properties. 

Injection  of  an  underexpanded  sonic  or  supersonic  jet  into  a  supersonic  crossflow 
produces  several  flow  structures.  Figure  2.1  shows  some  of  these  features  described  below. 
First,  a  bow  shock  is  produced  as  the  freestream  impacts  the  injection  streamtube  as  if  the 
injectant  were  a  cyhndrical  body.^  A  separation  bubble  can  form  slighdy  upstream  of  the 
injector  port  when  7  is  on  the  order  of  one  or  more.  This  is  due  to  the  boundary  layer 
being  thicker  than  the  injector  diameter.  After  entering  the  freestream,  the  jet  experiences  a 
rapid  Prandtl-Meyer  expansion  surrounded  by  a  barrel  shock.  A  shock  that  is  normal  to 
the  jet  path  terminates  the  barrel  shock  and  recompresses  the  flow.  This  shock  is  known  as 
the  Mach  disk  and  large  scale  mixing  occurs  after  this  normal  shock. 
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Figure  2.1:  Example  injection  schlieren 


2.2  Transverse  Injection 

Transverse  injection  is  the  simplest  injection  configuration.  It  requires  only  an  injection 
port  at  the  desired  injection  angle,  usually  normal  to  the  plane  of  the  injection  wall  and 
parallel  with  the  direction  of  flow.  Figure  2.2  illustrates  the  simple  transverse  injection  used 
in  this  study. 


Figure  2.2:  Simple  transverse  injection 
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The  entry  of  the  injectant  is  treated  as  a  two  stage  process.*'  The  jet  remains  relatively  intact 
as  it  first  enters  the  freestream  and  then  expands  to  the  height  of  the  Mach  disk.  Beyond  the 
Mach  disk,  the  flow  turns  and  accelerates  with  the  main  flow.  In  the  second  stage  the  jet 
acts  as  a  coaxial  vortex  mixing  structure. 

The  dynamic  pressure  ratio,  q  ,  strongly  impacts  mixing  and  penetration  in  the  near¬ 
field.  The  initial  rate  of  mixing  is  proportional  to  —  This  can  be  interpreted  from 

continuity  and  momentum  principles.  As  the  dynamic  pressure  of  the  injectant  is  increased, 
the  injectant  will  initially  enter  the  freestream  either  faster  or  denser.  If  the  injectant  were 
faster,  the  jet  itself  would  reduce  the  residence  time  it  has  to  mix  in  the  local  freestream.  If  it 
were  denser,  it  would  require  greater  time  and  distance  to  disperse.  Penetration  is  enhanced 
with  increasing  q  in  that  the  Mach  disk  is  pushed  further  out  into  the  freestream.^ 

2.3  Vorticity  Generation 

Vorticity  is  the  driving  mixing  mechanism  in  the  near  field.^  Streamwise  vorticity  is 
generated  by  two  mechanisms  in  transverse  injection.  The  first  is  the  external  pressure 
differential  created  across  the  front  and  tear  faces  of  the  emerging  jet.  The  main  flow  must 
split  around  the  injectant  jet,  creating  high  pressure  regions  on  the  front  and  sides  of  the  jet. 
The  resulting  low  pressure  on  the  rear  of  the  jet  becomes  an  entrainment  region  for  the  local 
freestream.  This  results  in  a  counter-rotating  vortex  pair  inside  the  injectant  jet.  The  second 
mechanism  is  created  as  the  plume  is  bent  coaxially  with  the  freestream.  The  inside  region 
of  the  plume  which  is  not  as  greatly  accelerated  is  drawn  into  the  center  of  the  plume. 

For  injection  over  a  ramp,  the  primary  vortex  generation  comes  from  the  ramp,  not  the 
injectant  plume.  Flow  pressurization  occurs  in  the  shock  layer  between  the  ramp  induced 
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bow  shock  and  the  ramp  compression  face.  The  high  pressure  fluid  then  spills  around  the 
ramp  causing  a  streamwise  vorticity.^^  The  vorticity  causes  large  scale  convective  mixing 
that  forms  a  counter-rotating  vortex  pair.^  The  relative  strength  of  the  vorticity  generated 
depends  on  the  freestream  Mach  number.  Additional  vorticity  is  imparted  by  the 
pressuri2:ed  flow  expanding  over  the  rear  face  of  the  ramp  and  wrapping  around  the  low 
pressme  region  encircling  the  underexpanded  jet  plume. 

Once  the  counter-rotating  vortex  pair  is  generated,  it  is  the  dominant  mixing  mechanism 
in  the  near  field.^’®  A  cross  section  view  of  the  developed  plume  remains  circular  with  the 
vortex  structures  clearly  visible  with  intense  rotation  rates.^^ 

2.4  Gas  Sampling 

For  the  complete  description  of  a  flowfield  generated  by  a  fully-constrained  supersonic 
flow  with  foreign  gas  injection,  it  is  necessary  to  obtain  species  composition  measurements. 
Gas  injection  creates  a  mixture  of  different  species  with  different  ratios  of  specific  heats,  y. 
Nearly  all  of  the  equations  relating  thermodynamic  variables  to  Mach  number  involve  y  as 
an  exponential  power.  Therefore,  only  a  slight  variation  may  cause  significant  changes  in 
reduced  thermodynamic  values. 

The  major  addition  to  this  study  over  previous  studies  was  the  different  gas  injectants 
used.  Because  of  this,  the  method  of  finding  the  concentration  levels  and  therefore  y  will  be 
discussed  in  detail  here.  The  concentration  measurement  probe  was  designed  to  include  a 
hot-fihn  sensor,  pressure  tap,  and  thermocouple.  The  condition  for  mass  continuity  is 
expressed  as 
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Since  the  flow  through  the  analyzer  is  choked  (M=l)23,  this  continiaity  equation  reduces  to 

P.C* 


where  C*  is  the  critical  flow  function  defined  by 
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where  the  *  designates  sonic  conditions  at  the  orifice  and  3  is  the  molecular  weight.  Since 
the  mass  flow  through  the  sensor  plane  must  equal  the  mass  flow  through  the  choked 
orifice,  the  mass  flux  at  the  sensor  plane  may  be  expressed  as 

P,C*  A* 


pu  = 


(2.4) 


Hence,  the  mass  flux  through  the  sensor  plane  is  a  function  of  the  total  pressure,  total 
temperature,  area  ratio,  and  gas  composition. 

Now  consider  the  heat  transfer  from  the  hot-film  sensor  at  the  measurement  plane  to 
the  passing  fluid.  The  rate  of  heat  transfer  is  given  by 

q,  =  ijR,  (2.5) 
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where  R^is  the  filra  resistance  and  J^is  the  film  current  supplied  by  the  anemometer. For  a 
constant  temperature  anemometer,  the  film  current  is  related  to  the  anemometer  response 
voltage  by 


V 

Rf+Rs 
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where  R^is  the  series  bridge  resistance.  Now  define  the  Nusselt  number  in  terms  of  the  rate 
of  heat  transfer  from  the  hot-film  and  the  temperature  difference  between  the  film  and  the 
passing  fluid  as 


Nu^ 


(2.7) 


where  k  is  the  thermal  conductivity  of  the  gas  mixture  and  /  is  the  active  sensing  length  of 
the  hot-film.23  Then,  combining  Equations  2.4,  2.5,  and  2.6,  the  Nusselt  number  may  be 
expressed  as 


Nu  = 


R. 
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(2.8) 


The  Nusselt  number  may  also  be  related  to  the  mass  flux  via  the  Reynolds  number. 
Re  =  ^ ,  where  d  is  the  film  diameter,  as 
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and  a  and  b  are  empirically  determined  constants.  Thus  equating  equations  2.7,  2.8,  and  2.9 
for  Nusselt  number  and  inserting  Equation  2.4  for  mass  flux  at  the  sensor  plane,  the 
governing  equation  for  the  anemometer  response  voltage  is  obtained  as 


(2.10) 


The  parameters  C*,  and  //  are  properties  of  the  gas  mixture  and  are  therefore  a 
function  of  the  gas  concentration,  Furthermore,  since  a  constant  temperature 

anemometer  is  employed,  the  film  resistance  T)and  therefore  R^are  constant.  Finally,  the 
values  of  /,  and  are  system  parameters  and  remain  constant.  Therefore  the  anemometer 
response  voltage  is  a  function  of  the  total  pressure  ,  the  total  temperature  T),  and  the  gas 
concentration,  That  is 

V=f(p,T,X^)  (2.11) 

Thus  for  a  given  total  pressure  and  total  temperature,  the  gas  concentration  is  uniquely 
determined  by  the  anemometer  response  voltage.  It  is  this  voltage  that  is  recorded  for  each 
concentration  and  pressure  dxaring  the  caUbration  and  data  runs.  FuUer  presents  a  detailed 
derivation  of  the  concepts  presented  in  this  section.^^ 
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3.  Facilities  and  Instrumentation 


This  chapter  contains  descriptions  of  the  facihties  and  instrumentation  used  during  this 
investigation.  Descriptions  of  the  AFIT  Mach  three  wind  tunnel,  the  transverse  injector, 
PME  ramp,  mean  flow  probes  and  gas-sampUng  instrumentation  ate  presented  followed  by 
descriptions  of  the  data  acquisition  system. 

3.1  AFIT  Mach  Three  Wind  Tunnel 

The  tests  for  this  investigation  were  conducted  using  the  AFIT  Mach  three  wind  tunnel. 
This  tunnel  was  a  blow-down  design,  using  a  pressure-vacuum  system.  The  budding- 
suppUed  pressurked  air  was  provided  by  two  Adas  Copco  GAU  807  compressors  at  a  0.5 
kg/s  flow  rate.  The  air  was  dried  by  two  Pioneer  R500A  refrigerant  air  dryers  en  route  to 
the  pressure  side  of  the  tunnel.  Final  drying  and  filtering  was  performed  by  a  cyclone 
separator  and  multiple  layers  of  Fdtrite®  filter  paper. 

Sixteen  vacuum  tanks  of  16  m^  volume  each  comprised  the  vacuum  system.  The 
evacuation  was  accomphshed  by  three  Stokes  Micro  Vac  pumps  and  the  tanks  were 
evacuated  to  approximately  10.0  mm  Hg  prior  to  each  run.  This  provided  over  twice  as 
much  run  time  as  was  needed  during  this  investigation. 

The  tunnel  itself  consisted  of  a  plenum  chamber  upstream  of  a  converging-diverging 
noz2le  to  produce  supersonic  flow.  Flow  straightening  screens  were  placed  between  the  air 
supply  and  plenum  chamber.  An  Endevco  0-960  kPa  pressure  transducer  mounted  in  the 
plenxom  chamber  measured  total  gauge  pressure.  Fluctuation  in  total  pressure  during  test 
runs  were  no  more  than  ±1.5%.  An  Omega  Engineering  type  K  thermocouple  provided 
plenum  total  temperature.  Fluctuations  in  total  temperature  during  test  runs  were  no  more 
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than  ±0.5%.  The  nominal  Mach  number  of  the  tunnel  no22le  is  2.8,  ±1.8%  with  a  standard 


operating  total  temperature  of 

3.2  Transverse  Injector  Model 


The  injector  model  is  described  in  full  detail  by  McCann.^o  The  injector  setup  was 
originally  used  to  research  turbulence  at  low  angles  of  injection.  Wilson  made  minor 
modifications  to  the  injector  model,  incorporating  PME  ramps  to  study  the  enhancements 
of  mixing  associated  with  the  ramps  and  to  provide  smoother  flow  from  the  injector.^^ 

Since  the  primary  objective  of  this  study  is  to  measure  differences  in  penetration  and  mixing 
of  different  injectant  gasses,  no  further  modifications  were  made  to  the  transverse  injector. 
Figure  3.1  is  a  side  view  of  the  test  section. 


Direction  of  Flow 


► 


Probe  measurement 
location,  x/d=20 


Pressurized  Injectant 

Figure  3.1:  Side  view  of  test  section 
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The  throat  diameter  of  the  injector  outlet  was  0.32  cm  expanding  to  0.40  cm  across  the 
exit  minor  axis.  The  settling  chamber  was  monitored  by  an  Endevco  0-690  kPa  pressure 
transducer.  A  high  pressure  line  was  attached  through  a  regulator  from  the  gas  bottles  to 
supply  argon  and  hehum.  The  gas  first  traveled  through  2  meters  of  copper  tubing  to  ensure 
the  injectant  gas  was  at  the  same  total  temperature  as  the  freestream.  An  Omega 
Engineering  type  K  thermocouple  was  inserted  in  the  settling  chamber  of  the  injector  to 
monitor  this  temperature.  Table  3.1  summarizes  the  injector  jet  conditions.  The  variations 
in  the  following  parameters  were  no  more  than  2.3%. 


Table  3.1.  Injector  jet  flow  conditions 


Mach 

1.82 

Poi 

304541  Pa 

Toi 

294  K 

Pi 

47595  Pa 

Ti 

173  K 

U.(He) 

1409.8  m/s 

Ui  (Ar) 

445.9  m/ s 

Y(He) 

1.67 

*  Y  (Ar) 

1.67 

>.(He) 

1.24 

?t(Ar) 

3.93 

9  (He) 

2.89 

^(Ar) 

2.89 

It  is  interesting  to  note  that  q  is  the  same  for  both  hehum  and  argon.  This  is  due  to  the 
differences  in  density  and  velocity  for  each  of  these  gases.  While  hehum  has  a  low  density, 
the  velocity  is  high,  which  makes  the  dynamic  pressure  ratio  with  constant  injector  pressure 
equal  to  argon  with  a  high  density  but  low  velocity. 
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3.2.1  Tunnel  Coordinate  System  and  Dimensions 


The  tunnel  coordinate  system  relative  to  the  injector  nozzle  exit  is  summarized  in  Figure 


3.2. 


Figure  3.2:  Tunnel  coordinate  system 

The  physical  dimensions  of  the  tunnel  and  injector  are  summarized  in  Table  3.2  below. 


Table  3.2.  Wind  tunnel  and  injector  nozzle  exit  dimensions 


Test  Section  Width 

6.35  cm 

Test  Section  Height 

6.35  cm 

Test  Section  Length 

35.0  cm 

Nozzle  Exit  Minor  Axis  Diameter,  d 

0.40  cm 

Nozzle  Exit  Major  Exit  Diameter,  dm 

0.92  cm 

Distance  from  Tunnel  Nozzle  to  Test 
Section 

27.0  cm 

Position  of  Injector  Exit  on  Test  Section 

8.5  cm 

(from  test  section  leading 
_ Edge) _ 

The  standard  freestream  conditions  for  a  nominal  test  run  are  given  in  Table  3.2. 


Table  3.3.  Freestream  flow  conditions 


Mach 

2.9 

Po 

217875  Pa 

To 

294  K 

p 

8028  Pa 

T 

114K 

Uo 

608  m/ s 

3.3  PME  Ramp 

Along  with  a  flat  plate,  a  single  penetration  and  mixing  enhancement  ramp,  the 
symmetric  ramp,  SR  3,  used  by  Wilson,  was  chosen  for  the  study. 22  This  ramp  was 
symmetric  about  the  plane  of  the  tunnel  coordinate  system.  This  ramp  introduced  a  two- 
dimensional  shock  because  of  the  flat  plate  compression  surface.  The  plume  interacts  with 
this  shock  while  the  trailing  edge  accelerates  the  flow  behind  the  ramp.  See  Wilson  for 
complete  design  criteria  and  dimensions  of  this  ramp. 

3.4  Measurement  Locations 

The  measurement  location  for  the  mean  flow  data  was  at  x/ d=2Q.  Because  of 
symmetry,  the  ^  direction  horizontal  coverage  of  the  tunnel  consisted  of  10  stations 
spanning  half  of  the  width.  Each  station  was  0.138  cm  apart  for  an  effective  coverage  of  4.0 
cm.  In  the  tunnel  coordinate  system,  this  was  a  range  from  -4.05  to  4-4.05  ^ d.  This  allowed 
the  same  spatial  resolution  as  Wilson  for  mean  flow  measurements.^^  Vertically,  the  span 
covered  was  2.54  cm  or  aj// d  from  0.3  to  6.3. 

3.5  Mean  Flow  Measurements 

Mean  flow  measurements  were  acquired  through  a  system  consisting  of  pressure  probes, 
transducers  and  a  traverse  mount  described  in  the  following  sections. 

3.5.1  Pressure  Probes 

A  Pitot  and  cone-static  probe  were  used  to  collect  the  mean  flow  data.  The  Pitot  probe 
was  constructed  of  1.6  mm  outer  diameter  stainless  steel.  The  pressure  was  fed  by  Tygon 
tubing  to  an  Endevco  0-103  kPa  transducer.  The  cone-static  probe  consisted  of  four  evenly 
spaced  0.34  mm  ports  that  fed  into  a  common  chamber.  This  design  has  been  shown  to 
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allow  a  misalignment  error  of  up  to  6°  off  centerline  of  the  flow  without  an  effect  on  the 
measured  pressure .22  The  cone-static  probe  pressure  was  also  fed  through  Tygon  tubing 
into  an  Endevco  0-103  kPa  transducer.  The  output  of  the  transducer  for  both  probes  was 
filtered  by  an  Endevco  model  4225  signal  conditioner.  The  signal  conditioner  was  assigned 
and  calibrated  for  each  transducer  used. 

3.5.2  Probe  Traverse  System 

The  pitot  and  cone-static  probes  were  swept  through  the  vertically  orientedj// d  range  for 
each  of  the  10  ^axis  stations  across  the  tunnel  cross  section.  The  traverse  system  used  was  a 
Dantec  3-axis  traverse  system.  The  position  of  the  traverse  was  recorded  as  a  voltage  by  a 
TransTek  Model  0217  linear  voltage  displacement  transducer.  The  voltage-position  relation 
was  linear  for  the  range  of  motion  of  the  probes. 

3.5.3  Data  Acquisition 

Data  acquisition  for  the  flow  measurements  was  accomplished  through  a  Dantec  90N10 
based  data  acquisition  board  connected  to  an  Acer  Pentium  computer.  Dantec 
Streamware®  software  acquired  the  raw  voltage  readings  collected  by  the  pressure  probes. 
The  sample  rate  used  for  data  acquisition  on  all  channels  was  400  Hz.  The  data  were 
converted  to  tab  delimited  text  files  for  processing  by  the  data  reduction  routines. 

3.6  Gas  Sampling 

Gas  sampling  for  concentration  measurements  were  obtained  using  an  inttusive 
sampling  probe  and  gas  analyzer  designed  for  supersonic  flow  by  Fuller.^^ 
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3.6.1  Sampling  Probe 

The  probe  used  in  the  current  study  was  first  proposed  by  Ninnemann It  consisted  of 
a  small  Pitot  tube  with  conical  tip  and  an  internal  divergence.  This  ensured  an  isokinetic 
sample  by  swallowing  the  shock  into  the  probe.  The  probe  was  constructed  from  stainless 
steel  with  an  outer  diameter  of  3.18  mm.  The  sampling  probe  was  approximately  25  cm 
long  and  was  terminated  by  the  gas  analyzer.  The  probe  was  attached  to  the  same  traverse 
system  described  previously  and  measurements  were  taken  over  the  same  grid. 

3.6.2  Gas  Analyzer 

The  gas  analyzer  consisted  of  a  rectangular  (84mm  x  36  mm  x  19  mm)  stainless  steel 
block  housing  with  a  3.86  mm  bore  through  the  center.  The  sensor  plane  consisted  of  a  TSI 
1220-20  high  temperature  platinum  hot-film  sensor  mounted  perpendicular  to  the  flow 
entering  the  analyzer,  a  Genisco  PB-923  pressure  transducer  with  a  range  of  0  to  100  psia, 
and  an  Omega  Engineering  Type-K  thermocouple.  The  flow  then  exited  the  probe  through 
a  choked  orifice. 

3.6.3  Electronics 

Data  acquisition  was  performed  with  a  National  Instruments  AT-MIO-16X 
multifunction  input/output  board  installed  in  a  486  rack-mounted  computer.  The  AT-MIO- 
16D  has  a  lOpsec,  16-bit,  sampling  ADC  that  can  monitor  up  to  8  different  channels  with 
programmable  gains. 

The  hot-film  sensor  was  used  in  conjunction  with  a  Dantec  Type55M01  constant 
temperature  anemometer  fitted  with  a  DISA  Type-55M10  CTA  standard  bridge.  The  output 
from  the  sensor  via  the  anemometer  was  passed  through  a  Frequency  Devices  Model  9002 
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dual-channel  programmable  filter.  The  filtered  signal  was  input  into  the  computer  via  the 
AT-MIO-16X  I/O  board. 

The  pressure  excitation  voltage  and  output  gain  was  suppUed  by  a  Measurements  Group 
23120  signal  conditioning  amphfier.  The  voltage-to-pressure  cahbration  was  obtained  via  a 
Druck  DPI  500  pressure  indicator.  The  output  voltages  for  pressures  ranging  from  0  to  100 
psia  and  temperatures  ranging  from  482  R  to  1100  R  were  recorded  using  the  above 
mentioned  data  acquisition  system. 

3. 7  Flow  Visualization 

Color  schheren  photographs  were  taken  for  each  gas  and  ramp  configuration.  Due  to 
the  limited  coverage  provided  by  the  optical  glass  side  walls  of  the  wind  tunnel,  two  views 
were  needed  to  form  a  composite  image  of  the  flowfield  around  the  injector  and  ramp.  The 
Hght  source  used  was  a  Cordin  Model  5401  arc  Hght  and  ahgned  by  hand.  The  Ught  was 
collimated  by  a  153  cm  focal  length  mirror  and  passed  through  the  test  section.  A  Nikon 
N90s  camera  with  a  Kodak  DCS420  digital  attachment  was  used  to  capture  the  images.  The 
lens  was  a  Tokina  80-200  mm  telephoto  with  a  Vivatar  2x  macro  focusing  teleconverter. 

3. 8  Other  Computational  Resources 

Data  reduction  was  performed  on  an  Ameritest  Pentium  PC.  Contour  plots  were 
generated  on  the  AFIT  Sun  Sparc  2  and  10  workstations  using  Tecplot  7.0.  The  schheren 
images  were  generated  by  Adobe  Photoshop. 
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4.  Data  Reduction 


This  chapter  details  the  data  reduction  routines  applied  to  the  raw  pressure  probe  data 
collected  with  the  cone-static  probe,  the  pitot  probe,  and  the  concentration  probe. 

4.1  Mean  Pressure  Measurements 

For  both  flow  pressure  probes,  the  raw  data  were  converted  from  gauge  pressures  to 
absolute  pressures  and  reduced  to  mean  pressure  measurements.  This  was  accomplished 
using  a  combination  of  VISUAL  BASIC  data  reduction  routines  and  EXCEL  spreadsheets. 
The  routines  first  scaled  the  probe  voltages  to  psi  by  the  appropriate  factor  depending  on 
probe  type.  Second,  the  gauge  psi  values  were  converted  to  inches  of  mercury  and  added  to 
the  atmospheric  pressure  in  the  wind  tunnel  lab  at  the  time  of  the  data  run.  Third,  the 
converted  pressures  were  averaged  by  groups  of  44  points  into  a  single  point.  The 
spreadsheet  also  determined  the  position  of  the  probes  by  averaging  the  same  range  of 
traverse  position  travel  voltages  into  a  single  point.  At  the  400  Hz  sample  rate,  44  samples 
comprised  0.11  seconds  of  traverse  movement.  This  span  of  movement  covers  0.03  cm  for 
each  mean  datum  and  provides  a  spatial  resolution  of  80  mean  datum  over  the  2.54  cm  of 
traversal. 

4.1.1  Pressure  Ratio 

The  ratio  of  the  cone-static  and  pitot  pressure  data  (PJP,^  was  determined.  Because  of 
the  time  required  for  each  probe,  gas,  and  ramp  configuration,  the  cone-static  and  pitot 
pressure  data  were  collected  on  separate  wind  tunnel  runs,  and  on  separate  days  in  some 
instances.  To  eliminate  the  error  caused  by  measuring  the  pressures  on  different  runs,  the 
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probe  pressure  ratios  were  normalked  to  each  other  with  the  recorded  total  pressure  of  the 
flow  for  the  data  run  as  shown  in  equation  4.1. 

A  =  A  /:5l  (4.1) 

P  P  P 

•* /o  /  ■'  /o 

4.2  Gas  Concentration 

First,  a  calibration  of  the  gas  concentration  was  performed.  This  consisted  of  two  steps. 
A  static  ak-gas  mixture  of  known  concentration,  temperature  and  pressure  was  analyzed  at 
different  concentration  level  and  pressure  levels.  Then  the  data  was  reduced  using  a 
VISUAL  BASIC  code  to  find  the  constants  a  and  b  from  Equation  2.10.  The  calibration 
procedure  for  the  gas  analyzer  also  yielded  a  curve  fit  for  the  various  given  gas 
concentrations  and  a  fixed  calibration  temperature.  Figure  4.1  is  the  calibration  curve  for 


argon. 


Figure  4.1:  Argon  calibration  curves 
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These  ciirve  fits  gave  the  hot-film  voltage  as  a  function  of  sampling  ptessure.  Thus,  for  a 
pair  of  values,  hot-film  voltage  and  sampling  pressure,  the  gas  concentration  could  be 
interpolated  provided  the  sampling  temperature  was  close  to  the  calibration  temperature. 

For  example,  say  that  at  a  given  sampling  pressme  and  temperature,  an  anemometer 
response  voltage  of  F  is  measured.  On  the  calibration  map,  at  the  measured  sampling 
pressure,  this  voltage  is  bounded  by  a  higher  voltage,  V*,  representing  a  higher  gas  mole 
fraction  of  X*"  and  a  lower  voltage,  V',  representing  a  lower  gas  mole  fraction  of  X'.  So  that 
the  linear  interpolation^^  will  yield  the  measured  gas  mole  fraction  as 

X  =  (4.2) 


4.3  Mach  Number 


To  determine  the  local  Mach  number,  the  pressure  ratio  and  gas  mole  fraction  were  used 
as  follows.  First,  the  specific  gas  constant,  R,  for  the  local  gas  concentration  was  found  by 
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R  = 


(4.3) 


Next  a  ratio  of  specific  heats  was  determined  by 


(4.4) 


Y  = 


(4.5) 


Finally,  the  Mach  number  was  determined  from  the  ratio  of  cone-static  to  Pitot  pressure. 
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Py P^2  fusing  a  table  look  up.  For  supersonic  flow  over  a  cone,  the  governing  equation 
was  derived  by  Taylor.^^  This  equation  is  a  second  order,  non-linear,  ordinary  differential 
equation  which  requires  a  numerical  solution.  Taylor  showed  that  the  ratio  of  cone-static 
surface  pressure,  to  the  freestream  pressure,  P;,  is  a  function  of  the  freestream  Mach 
number,  and  y.  This  ratio  divided  by  the  Rayleigh-Pitot  formula  yields  the  Mach  number 
as  a  function  of  the  ratio  of  cone-static  to  Pitot  pressure,  PJ P^2  ^^d  y.  Thus,  for  a  given 
cone-static  pressure,  Pitot  pressure,  and  y,  the  corresponding  Mach  number  may  be 
bilinearly  interpolated  using  an  appropriate  look-up  table.  The  Taylor-McCoU  equation  was 
solved  for  a  large  range  of  Mach  numbers.  These  solutions  were  combined  with  solutions  of 
the  Rayleigh-Pitot  formula  to  form  the  table. 
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5.  Results  and  Analysis 

This  chapter  describes  the  flow  field  for  each  of  the  four  configurations  (two  gases  and 
two  ramp  configurations).  For  all  configurations  studied,  the  following  order  of  analysis  is 
used.  First,  schlieren  flow  visualization  images  are  used  for  qualitative  results.  The  basic 
flow  shock  structures  along  the  centerline  x-j  plane  are  examined.  Second,  the  gas  mole 
fraction  concentration,  Mach  number  contour  plots,  and  total  pressure  ratio  contour  plots 
are  presented.  The  contour  plots  are  presented  as  two  half  plots.  Since  the  flow  is 
symmetric  about  the  x  axis,  only  half  of  the  flow  is  needed  for  each  case.  This  allows  easy 
comparison  between  two  cases.  Finally,  the  apparent  effects  of  molecular  weight  are 
discussed. 

5.1  HeUum  Injection 

5.1.1  Flow  Visualization 

Figure  5.1  is  the  schlieren  photograph  of  helium  injection  with  no  ramp.  Figure  5.2 
shows  the  case  of  helium  injection  with  a  ramp  present.  All  of  the  shock  structures 
discussed  in  section  2.1  can  be  seen  in  both  photographs.  Note  the  shock  created  parallel  to 
the  ramp  face  in  figure  5.2,  Also,  the  expansion  of  the  flow  behind  the  ramp  is  apparent. 
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Figuae;  5.1 :  ScWieuen  photogtaph  of  helium  mjection  uath  no  lamp 


Figuae:  5.2:  Schlieaen  photogiaph  of  helium  injection  uath  lamp 
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5.1.2  Concentration 


Figure  5.3  presents  the  molar  concentration  data  for  helium.  The  with-ramp 
configuration  has  a  smaller  area  of  concentrated  helium  at  a  core  of  the  vortex  while  the  no¬ 
ramp  case  has  a  large  area  of  high  concentration.  The  ramp  helps  to  disperse  the  helium  to  a 
wider  area  of  the  plume,  while  the  no-ramp  case  concentrates  the  helium  within  the  plume. 
Also  the  plume  area  was  computed  as  17%  larger  with  the  ramp  than  without.  This 
compares  well  with  a  30%  increase  in  plume  size  for  injection  with  air.^ 
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Figure  5.3;  Helium  concentration  contour 


5.1.3  Mach 

As  expected,  the  Mach  contour  shown  in  Figure  5.4  shows  a  low  Mach  number  within 
the  core  of  the  plume  for  the  ramp  case,  while  without  the  ramp  the  plume  has  higher  and 


more  uniform  Mach  numbers.  The  plume  size  was  found  to  be  23%  larger  with  the  ramp. 


giving  consistent  results  with  the  plume  area  measured  with  concentration.  This  is  also 


300%  the  area  of  the  plume  found  with  air  injection. ^2 
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Figure  5.4:  Helium  Mach  contour 


Figure  5.5:  Helium  total  pressure  ratio 


5.1.4  Pressure 


Figure  5.5  shows  a  comparison  of  the  total  pressure  ratio  with  helium  injection  with  and 
without  the  PME  ramp.  The  ramp  clearly  increases  the  ratio  within  the  plume.  While  the 
pressure  loss  from  the  ramp  is  not  desirable,  the  additional  mixing  that  is  attained  because  of 
the  ramp  as  seen  in  the  Mach  and  concentration  contours,  outweighs  this  drawback.  This 
data  is  consistent  with  prior  study  results.^^  Compared  to  the  simple  traverse  injection,  the 
ramp  results  indicate  the  vorticity  of  the  plume  has  been  substantially  diminished,  but  mixing 
and  penetration  is  enhanced. 

5.2  Argon  Injection 

5.2.1  Flow  Visualization 

Figure  5.6  is  the  schlieren  photograph  of  argon  injection  with  no  ramp.  Figure  5.7 
shows  the  case  of  argon  injection  with  a  ramp  present.  Again,  aU  of  the  shock  structures 
discussed  in  section  2.1  can  be  seen  in  both  photographs.  Note  the  shock  created  parallel  to 
the  ramp  face  and  expansion  behind  the  ramp  seen  in  figure  5.2  for  helium  are  also  present 
in  figure  5.7. 
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Figuxe  5.6:  SchlieiBn  photogcaph  of  aigon  injection  \mth  no  camp 


Figuce  5.7:  Schlieten  photogcaph  of  aigon  injection  vdth  camp 
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5.2.2  Concentration 


Figure  5.8  shows  the  molar  concentration  data  for  argon.  Again,  we  see  that  the  no 
ramp  case  has  a  larger  area  of  high  concentration  argon.  Both  cases  concentrate  the  argon 
around  the  outer  edge  of  the  plume.  This  was  not  seen  with  helium  and  is  attributed  to  the 
higher  molecular  weight  of  argon,  where  the  argon  is  spun  to  the  outside  of  the  plume  by 
the  vortex.  The  overall  plume  area  for  the  with-ramp  case  is  26%  larger  than  without.  This 
is  very  similar  to  the  30%  seen  with  air.22  This  can  be  attributed  to  the  similar  molecular 
weight  for  argon  and  air.  Helium  is  much  lighter  than  air  or  argon  and  therefore  the  gas 
mixes  into  the  freestream  faster.  Air  and  argon  are  relatively  denser,  slower  moving  gases 
and  are  kept  within  the  plume  for  a  longer  time. 
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Figure  5.8:  Argon  concentration  contour 


5.2.3  Mach 

Again  as  expected,  the  Mach  contour  shown  in  Figure  5.9  shows  a  larger  area  of  low 
Mach  number  within  the  core  of  the  plume  for  the  ramp  case  while  higher  Mach  numbers 
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ate  seen  without  the  tamp.  The  plume  size  was  found  to  be  42%  latget  with  the  tamp, 
giving  consistent  tesults  with  the  plume  atea  measured  with  concenttation.  This  is  also 
about  the  same  size  plume  found  with  ait  injection.22  The  effect  of  moleculat  weight  can  be 
seen  in  figure  5.9.  Air  and  argon  have  similar  weights  and  the  Mach  contours  are  similar, 
while  the  lighter  helium  has  a  much  larger  plume  area  as  indicated  in  figure  5.4. 

5.2.4  Pressure  Data 

Figure  5.10  shows  a  comparison  of  the  total  pressrue  ratio  for  argon  injection  with  and 
without  the  PME  ramp.  The  ramp  clearly  increases  the  ratio  within  the  plume  as  was  also 
seen  with  hehum.  This  data  is  consistent  with  prior  study  results  as  expected.22 

5.3  Gas  Comparison 

In  both  the  ramp  and  no  ramp  cases,  hehum  provides  the  greatest  mixing.  The  Mach 
number  determined  plume  size  when  compared  with  air,  was  three  time  larger  with  hehum 
as  the  injectant.  Argon  produced  a  pliune  about  the  same  size  as  air.  This  larger  plume  size 
indicates  greater  mixing  and  penetration  of  the  plume  into  the  flowfield.  Therefore 
molecular  weight  does  effect  the  mixing  and  penetration  of  the  injectant  into  the  freestream. 
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6.  Conclusions  and  Recommendation 


This  chaptet  outlines  the  primary  conclusions  drawn  from  the  analysis  in  Chapter  Five. 
Recommendations  for  further  study  are  submitted. 

6.1  Conclusions 

As  previously  stated,  the  primary  objective  of  the  investigation  was  to  determine  the 
differences  molecular  weight  would  cause  on  mixing  and  penetration  of  a  gas  injecting  into  a 
supersonic  freestream.  The  experimental  results  show  a  trend:  decreasing  molecular  weight 
increases  mixing  and  penetration  in  a  simple  transverse  injection  with  a  PME  ramp. 

In  general,  the  concentration  and  Mach  plots  show  agreement  with  expected  results. 

The  PME  ramp  performed  as  had  been  demonstrated  in  prior  studies.  The  differences  in 
the  injectants  are  clear.  The  lower  molecular  weight  allows  for  greater  mixing  and 
penetration  of  the  plume  into  the  freestream. 

6.2  Recommendation 

Results  have  been  established  which  clearly  indicate  the  trends  associated  with  the 
effects  of  molecular  weight  on  mixing  and  penetration.  The  next  step  in  conducting  the 
experiment  is  to  continue  using  hehum  and  to  further  explore  the  use  of  PME  ramps. 
Hydrogen  is  the  most  Ukely  candidate  for  SCRAMjet  applications  due  to  the  cooling 
possibilities.  Therefore,  fully  exploring  the  effects  of  low  molecular  weight  would  be 
beneficial. 

Additional  methods  can  be  added  to  this  investigation  to  provide  a  more  complete  view 
of  the  flow  structures.  Laser  imagery  can  be  used  to  enhance  the  flow  visualization  to  better 
interpret  the  flow  data.  Also,  hot  wife  anemometry  can  be  used  to  direcdy  measure  the 
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strength  of  vorticity  and  turbulent  properties.  This  would  allow  complete  characterization  of 
the  flow. 

Finally,  a  continuation  of  the  techniques  employed  in  this  investigation  is  recommended. 
A  study  of  the  trends  of  molecular  weight  and  PME  effects  have  been  tested  and  evaluations 
made.  Future  research  can  be  focused  on  broadening  the  basic  understanding  of  the 
complex  flow  fields,  perhaps  in  the  presence  of  combustion.  Once  a  suitable  base  of 
knowledge  exists,  refinement  of  the  design  choices  of  fuel  and  mixing  enhancers  can  lead  to 
optimum  configurations. 
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Appendix  A:  Uncertainty  Analysis 


A.1  Position  Uncertainty 

The  measixrements  conducted  in  this  experiment  were  made  in  two  dimensions.  In  the 
horizontal  direction,  the  positions  for  measurement  were  determined  using  a  steel  ruler  and  a 
point  relative  to  the  tunnel  floor.  The  ruler  was  graduated  into  0.7938  cm  divisions.  A 
maximum  misalignment  of  Vi  of  the  axial  station  separation  distance  would  be  possible 
based  upon  the  ruler  graduation.  This  translates  into  a  maximum  z  axis  position  uncertainty 
of  ±0.08  cm. 

For  each  z  axis  station,  a  probe  sweep  of  2.45  cm  of  the  tunnel  cross  section  was  made 
by  the  traverse  system.  The  vertical  position  of  the  probe  was  tracked  by  the  LVDT  system 
described  in  section  3.5.2.  The  position  indicated  by  the  LVDT  voltage  had  a  recorded 
variance  of  ±0.13%.  This  error  is  deemed  insignificant  when  compared  with  the  horizontal 
position  error. 

A2  Pressure  Data  Uncertainty 

The  acquired  mean  flow  pressure  data  have  two  associated  uncertainties  related  to  gain 
uncertainty  and  gain  stability  across  the  measurement  range.  The  gain  uncertainty  was  0.5% 
and  the  stability  uncertainty  was  0.2%.  Including  digital  conversion  effects,  Huffman  found 
the  overall  pressure  uncertainty  was  0.54%.  McCann  found  the  Mach  number  error  for 
similar  measurements  to  be  1%P 
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A,3  Gas  Analyzer  Uncertainty 


The  uncertainly  of  the  concentration  measurement  was  previously  determined  by  Fuller 
to  be  ±0.2%  in  mole  fraction,  and  an  estimated  three  percent  uncertainty  in  the  calibration 
curve  fir  and  interpolation  routine.^^  This  results  in  a  total  uncertainty  of  ±0.05  in  the  mole 
fraction. 
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